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Presented herein are the results of a flight test program 
conducted to determine two of the longitudinal stability deriva- 
tives of the Cessna 140 airplane. These derivatives are the 


elevator power, C, $? and the demping in pitch, Omae° The flight 


test techniques employed to obtain these derivatives are outlined. 


The values of C,, 6 and Cnag thus obtained are compared with the 
values of these parameters computed from the theoretical formulae 
as a means of evaluating the steady state flight test procedures 
employed. Sources of error which arise in the experimental de- 
termination of the derivatives are analyzed and techniques of 


Minimizing these errors are indicated. 











OBJuC TIVE 


The objective of this investigation has been the determination 
of two longitudinal stability derivatives -- the elevator power (C,, s ) 
end the damping in pitch (Snag) -- of a typical light airplans, the 
Casena 140. By conparing the exnerimental valuss of these ceriva- 
tives dorived from flight tert deta with the values calculated by 
means of the theoretice] arproach, it is hoped thst some evaluation 
or the flight test techniques employed to cbtain the derivatives can 
be mede. The results may be of some value to the lignt airplane 


desimner in determining the required elevator power for airnlancs 


with high demring in pitch. 


The investigation was conducted under the direction of the 
Princeton University Aeronautical Engineering Department at Princeton 


Airport, Princeton, New Jersey, from January to lay, 1950. 








INTRODUCTION 


In recent years the increasing speeds and performance of aircraft 
heve necessarily focused a attention to the stability and control 
problems involved in securing desireble hendling qualities of the airplane 
in high speed flight. jn the past, much of the necessary stability and 
control data could be successfully predicted with good accuracy from wind 
tunnel tests. Today, however, the high subsonic, transonic, and super- 
sonic speeds under consideration definitely limit the usefulness of the 
wind tunnel in this field. The difficulties of design and construction 
and the prohibitive costs of such tunnels, coupled with the inaccuracies 
of small seale models, force the aerodynamicist to turn to other meens 
‘to obtain stability and control data. The other possibility is the air- 
plane itself and this hes, in fact, been the modern approach to the 
problem. 

In general, it is necessary to deduce the sirplane's aerodynamic 
charecteristics from the response to its controls. ‘ich study is now being 
given to the so-called frequency response of the airplane as a means of 
determining the stability parameters. In this approach the airplane is 
oscillated by an oscillatory input to a particular control and its steady 
state frequency response is recorded. In some cases the frequency response 
may be deduced from the transient response to a pulsed input to the control. 
The eerodynanic parameters are then calculated from the frequency response 
data. While these methods will undoubtedly be of increasing usefulness in 
the future, they are at present somewhat limited by the need for precise 
instrumentation and accurate data reduction techniques. Exact knowledge of 
the airplane's mass and inertia characteristics is required and this in 


itself is no small problen. 
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Whore the ®liegnt duration of the eirplano is not limited severely, 
os in the ease of practically ali. oresont dey cubsonic aireraft, certain 
Sicht testing teclniques for obtainine some of the stability parameters 
have been developed which are much simpler in theory and execution than 
the frequency response methods. In many cases these teehniquos ara not 
onlv more accurete but their value is enhanced by lower ecst, cimpler 
instrumentation, and more precise data reduction. It is with some of 
these procedures that this paper will doal. 

The techniques under consideration involve the deduction of the air- 
plane's elevator nower and its damning in pitch from the eo-celled trin 
curves taken in unaccelerated and accelereted flight. As a means of 
evaluating the methods, flight tests of this type were conducted on a 
Cesena 140 airplane, and these two stability derivatives were computed 
from the results. Comperison with results obtained in computing the de- 
rivatives from the theorstical equations shoulec yield some insight into 


the accuracy that may be expected in determining these derivetives from 


flight tests in this manner, 
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DESCRIPTION OF THE AIRPLANE 


The Cessna 140 is a single engine, high wing, _ place, 
personal type monoplane with external bracing and fixed, conventional 
landing gear. The wing is rectangular with rounded tips. It has be 
the normal configuration, single vertical tail, Frieze type ailerons, 
and is equipped with trailing edge plain flaps. 

The airplane is of semi-monocoque all metal construction exce pt 
the wings which are fabric covered. All control surfaces are of 
metal. Figure IV-1 of the appendix illustrates the appearance of 
the airplane as instrumented for these flight tests. 

The following general specifications and dimensions are taken 
fron the drawings and reports of the manufacturer except in the cese 
where they have been modified by the installation of special equip- 
ment. In these cases both the old and the new dimension or specifi- 


cation is given. 


ns a a - _— eee ee —— 
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firplane, general 

Manufacturer 

Type 

Recommended gross weight 

Center of gravity range: 
forward limit 
aft linit 

Overall length 

Height 

Maxirum allowable maneuvering load factor 
gross was eht 1460 lbs. 
flaps down 40° 

Wing 

Airfoil section 

S pan 

eon (total) 


Area (less ailerons) 


Aspect retio 
Taper ratio 
Chord 


Mean aerodynamic chord 
Length 


Distance of leading edge back of nose 
reference dzetum line 


Incidence 


Dihedral 
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Cessna Aircraft, Co. 
140 
1450 lbe. 


22.85 mac. 
30.0% mac. 
2506 20 ig 


Thea hae 


this 57 to “Le 26 


fl 097 to mwe 26 


394 ine 
159.29 sa. ft. 


145.21 sq 


6.76 

1.0 

60.5 in. 

59.02 in. 
56.58 in. 
flee 


Are 








~ 


Aileron 

Type 

Area 

Span 

Chord 

. Travel 
up (from noutral) 
down (From neutral) 

ning flaps 

1ype 

Area 

Peri 

Chora 


Trevel (down) 


korizontal.tail curface 


Airfoil section 


Area (including elevators) 


» pan 
Maximam chord 
Incidence 


Dihedral 


Elevator area (total, including tab) 


Elevator span 


Elevator travel 


up (from streamline with ctebilizer) 


dom (from streamline with stabilizer) 
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“odified Frieze 
14.00 eq. FE. 
74 in. 

14 in. 


© ie te TU. 


o.0 ane 

L0° 

macs COC9 
24.35 sq. ft. 
1Co"4n. 
41.4 in. 
2a 

0 

9.66 sq. ft. 
106 in. 

20" 

20° 











ee tab area 
Elevator trim tab’span 
Elevator trim tab mean chord 
uwlevator trim tab travel 
uy (fron elevator trailing edge) 
dowa (from elevator trailing edge) 
Vertical tail surface 
frea 
Fin area 
Span (to fuselage center line) 
nudder Lrea 
Fudder span (maxim) 
hudder trevel 
right (from streamline with fin} 
left (from etreamline with fin) 
Fuselage 
Maximum width 
Maximum height 
Length (tip of noso to tip of tail) 
Dngine (Continental) 
Type 
Number cf Cylinders 
Power Settings Condition 
Teke off 


nrte 
Cruise 
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0.695 SQ. fte 
36 in. 
Seed iy 


12.420 sq. ft. 
6.663 OG e Stic 
52e2 ine 


Oe oe Se ft. 


40.0 in. 
57.0 im. 


256.50 in. 


C 85 
4 


25°75 
2575 
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Propeller 
Manufacturer 
Type 


Diameter 
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censenich Bros. 
Vood, fixed vitch 


79 ine 
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INSTRUMENTATION 


The instrumentetion required to obtain the necessary data for this 
investigation consisted of an airspeed indicator, an altimeter, an accelero- 
meter, and an elevator position indicator. The individual instruments 
may be described as fol lOws': 

1. Airspeed Indicator: 

The airspeed was measured with a standard sensitive type airspeed 
indicator connected to a full swiveling pitot-static head. This instru- 
ment was calibrated by means of the speed course method and a calibration 
curve is contained in the appendix. 

2. Altimeter: 
The altitude was measured with a standard sensitive type al ti eevee 
3. Accelerometer: 

Since the range of accelerations encountered in the tests was small 
and since it was desired to measure increments of acceleration as small 
as 0.25 g, the scale of a standard accelerometer was considered to be 
inadequate. A special accelerometer was constructed and calibrated for 
this investigation. 

This accelerometer consisted of a glass tube approximetely 20 inches - 
long, a small coiled spring, a steel weight, and an attachment for con- 
necting the spring inside one end of the tube. Two edditional weights 
were constructed for use in calibration. These two weights were in the 
ratios of two and three to the weight to be used in the instrument. 

With the spring mounted in one end of the glass tube, the instrument 


wes calibrated by securing the tube in a vertical position, hanging each 








Le 
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of the three weights on the spring in turn and carefully marking 
on the tube the equilibrium position of the weight. In this manner 
the tube was calibrated for accelerations of one, two, and three 
"etst, three "g's" requiring an extension of the spring of about 
ten inches. Assuming the spring constant to be linear, the fractional 
"et positions were easily located and inscribed on the tube. 

A photograph of this accelerometer as mounted in the airplane 


may be seen in the Appendix. 


Elevator Position Indicator: 


The elevator position was measured by a 26 volt, 400 cycle a-c 
autosyn type transmitter. This transmitter was linked to the elevator 
control cable as close to the elevator horn as practicable. A cali- 
bration of the autosyn in terms of elevator angle is contained in 
the Appendix. 


x 


All data from the flight tests was visually read and manually recorded 


by the pilot and co-pilot. 
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WEIGHT AND BALANCE 


All flight tests were made for three center of gravity locations, 
namely, 24.6%, 27.6%, and 30.6% mac. The normal center of gravity 
location for the airplane less pilots is 27.4% for a gross weight of 
1100 pounds. The weight of the two pilots and parachutes was 335 
pounds located at the center of gravity in all three cases. The normal 
operating gross weight for this airplane is 1450 pounds. The airplane 
was ballasted to the desired centers of gravity by means of lead 
weights attached either to the tail wheel arm or to the engine mounts 
as required. As all tests were flown with the fuel tanks approximately 
full, and since all tanks are located near the normal center of gravity, 
no corrections to center of gravity position were made to account for 


fuel consumption. 
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RESULTS AND DISCUSSION 


J. GENERAL 


The flight test date required for computing the two desired stability 
derivatives, the elevator power and the damping in pitch, is essentially 
the same in both cases. Determination of the elevator power, Ons » requires 
curves of 6, versus V; at 1 g for several center of gravity locations. De- 
termination of the damping in pitch, Caos requires these same curves for 
at least one center of gravity location (preferably several) and, in addition, 
the curves must be obtained not only at 1 g, but for several values of 
acceleration, i.e., several values of n. 

The data was taken for three center of gravity locations as indicated 
under "Weight and Balance". The speed range over which the tests were run 
extended from 60 mph, Vy, to 120 mph, Vy. For the 1 g curves, deta was taken 
in steady unaccelerated flight at iHeratients in airspeed of 5 mph. For 
values of n greater than 1.0, the data was taken in accelerated flight in 
steady banked turns at increments of 10 mph in airspeed. 

The airspeed instrument error was negligible; therefore V; was equal to 
Vee All flight tests were made at a pressure altitude of 5000 feet and 
since the speed range was low, the effects of compressibility were considered 
to be negligible and V, was taken equal to V,. The flight test data is 
plotted in Figures 1, 2, and 3 as curves of dg versus V,. 

In order to eliminate the effects of the elevator tab on required 
elevator deflection, and also to eliminate tab corrections in the calculations, 


the elevator tab was trimmed to zero deflection throughout the tests. 
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The greatest source of error (as will be shown later) in determining 
the desired derivatives, came in through the variation in the power effects. 
Since the Cessna 140 is equipped with a fixed pitch propeller it was im- 
possible to maintain a constant power setting while varying the speeq. At 
the higher end of the speed range (100 - 120 mph) it was necessary to 
reduce power rather drastically to avoid over-speeding the engine. The 
RPH was held constant throughout the tests. The cruising RPM of 2400 was 
chosen as the setting and was held at this value for all runs. 

The succeeding sections treat the two derivatives under consideration 
in some detail in attempting to show how the value of the derivative is 
determined from reduction of the flight test data end how this value compares 
in magnitude with the theoretically computed derivetive. Sources of error 
are analyzed and finally, an evaluation of the technique for obtaining the 


particular derivative is given. 


II. THE ELEVATOR POWER, Cm 


ihe Theory: 
When an airplane is in equilibrium the sum of the moments acting on 


the airplane about any chosen reference line must be zero. If, for instance, 
the root leading edge of the wing is selected as the datum, or reference 
line, then the moment of the airplane's weight, W, about this datum is the 
product of the weight and the distance from the reference line to the center 
of gravity of the airplane (assuming the weight of the airplane to be con- 


centrated at the center of gravity). This moment, then, is given by: 


M ref a Xe.g, (1 
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Now, if a new center of gravity position is obtained by shifting 
some portion of the airplane's weight (ballast weight s Wy) to a position aft 
of the original center of gravity by a distance, d, the moment of the air- 
plane's weight about the reference line is given by: 
: 
| ~ 4 
Miner, YV Cg. > (2) 


Therefore : 


W - Soe = (vw -Wp) e.g, + Wy, (d + Xe.g, ) 
WW Kog., = WW. Keg. — Wh Xcq, +Whpid + Wy: X 6.9, 


W (xeo.-Xe9,) = Weed | (3) 


Obviously, then, the product of the airplane's weight and the change in 


center of gravity position must equal the applied moment, i.e.: 


W (Keg, = Keg, ) = M applied | (4) 
In flight, any moment acting on the airplane in the longitudinal plane 
must be balanced by a moment produced by the deflection of the elevator, ° 
That is, the sum of the applied moment and the moment due to elevator de- 
flection mst be equal to zero. Therefore, the magnitude of the pitching 
moment coefficient produced by the elevator per degree deflection is a 


measure of the elevator power. The Soudnt Peechigge by the eleweilor is given 


by: 





Zoe = dMe : Ade (5) 
— A &e 
Therefore : 
W (xeg, - Xe9,,) + dMe. Ade = O _ (6) 





dSe 
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At equilibrium the lift must be equal to the weight (L » W). Making this 


substitution and putting the equation in coefficient form by dividing by 
qsc: 


Ci (X42 - 4e9-) = — Cm 5 ‘Ade (7) 


Solving for Cn § : 


C, (Xeg,— Xc.a)) | (8) 


me > ASe 


where X,_ Be is now given in per cent of the mean aerodynamic chord. 


C 


2. Determination of C,, 5! 

Examination of the expression for c. 5 developed in section 1 above, 
indieates that Cm § can be computed (as a function of airplane lift coefficient) 
from data obtained in simple flight test procedures. This technique may be 
outlined as follows: 

1. Curves of §, versus are obtained for at least two center of 

gravity locations of the airplane. 

2 For a given CL. determine from these curves the =e in elevator 

angle required by the given shift in center of gravity. 
| 3. Compute Cy» using equation (8). 

It should be noted that the curves of §, shown in Figure 4 are plotted 
versus the airplane's lift coefficient which is, in general, not equal to the 
lift coefficient of the wing alone. The simple relation derived above for 
determining the elevator power, On § » assumes that CL. 2 Cle That this is 
not necessarily true is easily recognized from the fact that the lift of the 
tail is a definite contribution to the total lift of the airplane and varies 


with elevator deflection. For CL. to remain a constant as the center of 
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gravity position is shifted requires that for a change in tail lift due 

to elevator deflection a corresponding increase or decrease in Cy, must 

occur. The change in wing angle of attack necessary to produce this lift 
change must in turn be compensated for by a slight change in elevator de- 
flection. It has been shown, however, that the errors introduced by 

assuming Cy, = ee for normal configurations are not large and probably 

are well within the accuracy of the flight tests. (Reference 3). For this 
reason it was considered justifiable to assume Cy, Ao Cr throughout this 
investigation and, therefore, Fh i= calculated using the simplified approach 


as outlined above. 


3. Results: 
A curve of C, g versus Cy was computed in the above manner as shown in 
Figure 5. Since three equally spaced center of gravity locations were used 


in the flight tests, the values of be used in computing C,,.were taken as 


mS 
the average change in elevator angle between the curves of §, versus C; for 
the three center of gravity locations. Theoretically, for equal incremental 
changes in center oP gravity position these curves will be equally spaced. 
That is, the change in elevator angle required to maintain a given Cy; will 
be a constant for even increments of center of gravity shift. Also, theo- 
retical curves of &@ versus C,, for various center of gravity locations will 
intersect at a common elevator deflection at C; = 0. It will be noted that 
the experimental curves of Figure 4 do not follow the theory exactly, 
particularly when C; is less than 0.4. The curves obtained for the aft and 


mid center of gravity locetions tend to follow the theory closely. However, 
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for lift coefficients below 0.4 the curve for the forward center of 
gravity exhibits a tendency to break down. For this reason the elevator 
power, Cr, 5? was computed only in the range of airplane lift coefficient 
from 0.40 to 1.20. 

The discrepancy found in the forward center of grevity curve was 
also exhibited in previous flight tests of the Cessna 140. (Reference 2). 
The decrease in required elevator angle at the higher speeds may be due 
to one or more of several possibilities. It has been shown (reference 3) 
that as the-center of gravity is moved further and further forward, the 
elevator increment per unit center of gravity change increases. However, 
the deviation of the curve is rather abrupt below Che =- 0.40 and it is 
believed that the discrepancy may be due in whole, or in part, to certain 
aero-elastic effects not presently accounted for. (Stabilizer deflection 


is among the possibilities.) 


4. Theoretical Calculation of C, ng? 
For an airplane flying in equilibrium (propeller off) the balance of 


moments (in coefficient form) is given by: 





Ome... i Omar, a Cu “2 v Cm sus, . = Vv df: (9) 


(Reference 1.) Differentieting this expression with respect to Se: 


= (Ae 





dC, = — (| Vim d X+ 


d Se AK I+ dée 
Ake ee elevater effectweness (10) 
de 


“ Omg = — a4' Vine Te | (11) 











ee ——s Oy a 
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This expression contains no propeller or power terms. This omission 
produces no error for the propeller wind-milling condition but will give 
results somewhat low for high power, low speed flight because of the 
increase in dynamic pressure at the tail due to the slipstrean. 

Calculation of the theoretical value of the elevator power according 


to this equation is as follows: 


@ 
me 8: Ay s abe = _ ACG eae =0ioee 
St IAL x 24.35 


This value for the tail aspect ratio is used to obtain the 
slope of the 1if% curve of the tail, a,, from Ref. 1, Fig. 5-5. 


a, = 0.048 per degree 


t = ei tale & lee = 0.228 


a SP 
~ aie 159.29 x 59.02 


<1 
<1 
t 
ale 
ce 
~ 


Be ht: Assume " = 0.90 


m= ©: e = 92:66. = 0.398 
e Sy 2h,» 35 


Using this value of the ratio of elevator area to horizontal 
tail area the theoretical value of SR is obtained from Ref. 1, 
Fig. 5-33. 


C2 = 0.59 


Therefore: Cg = —* (0.048 x 0.428 x 0.90 x 0.59) 


a 
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de Comparison of Theoretical and Experimental Values for C,, a: 

The theoretically determined value for the elevator power (C,, ce 
~~ 0.0109) is plotted in Figure 5 for purposes of comparison with the 
experimental curve. The theory assumes C,, § to be a constant based on the 
elimination of power effects. As pointed out above, this assumption leads 
to obvious inaccuracies in the low speed, high power condition due to the 
increase in the so-called tail efficiency, N° The tail efficiency is 
defined as the ratio of q at the tail to the free stream q. Ons varies 
directly as Ne and hence increases as C, increases. 

This variation of Cy 5 with C; is quite evident in the experimental 
Curve » Ong varying from — 0.0135 at C; = 0.4 to —0.0167 at C; = 1.2. 
Comparing these values with the theoretical value for this derivative, the 
experimental value is found to vary from 71.5% higher than theoretical for 
the low speed, high power condition (Cy, = 1.2) to 23.8% higher than theo- 
retical at cruising speed (C; = 0.4). 

As pointed out previously, however, with a fixed pitch propeller it 
was impossible to attain constant power settings over the speed range in- ° 
vestigated. In order to avoid over-speeding the engine, power was rather 
sharply reduced at the higher speeds. Such reduction, obviously induces 
a greater variation in n, than would normally be expected if the power could 
have been held constant. If this variation in the tail efficiency were 
accurately known and applied in making the theoretical calculation of Cy Es 
the discrepancy between theoretical and experimental values would tend to 
be reduced and probably would approach a constant as the lift coefficient 


increases. ‘+4 varies somewhat linearly with C,. Hence, applying this 
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correction to the tail efficiency would not affect the linearity of the 
theoretical curve, 

The experimental determination of Crs as performed under the as- 
sumption that oie = Oly also leads te values for the elevator power that 
are somewhat high. However, the use of this assumpotion yields an error 


that is usually within the accuracy of the flight tests. 


6. Eveluation of Experimental Techniouve for Obtaining Cm : 


The stick-fixed trim curves are among the most accurate data that may 
be obtained from flight testing an airplane. These curves are usually 
reduced to plots of elevator angle versus lift coefficient, Chaos for several 
center of gravity locations. As the pilot technique required in these tests 
is simple, the accuracy of the data depends only on how accurately the 
elevator angle and the airspeed can be measured (assuming, of course, that 
the gross weight and the center of gravity location are accurately known). 
Excellent instrumentation is available for measuring these parameters. 

The greatest errors involved in obtaining the elevator power with this 
technique come in through the power effects and the veriation in Ne For 
an airplane equipped with a controllable pitehn propeller, the power can at 
least be held constant throughout the speed range of the investigation. 
Nevertheless, it is extremely difficult to accurately predict the variation 
in Ne and herein lies the greatest source of error. The problem is further 
complicated by the fact that for various airplane configurations the hori- 
zontal tail may be positioned such that it is completely in, partially in, 


or conpletely out of the slipstream. In the latter case, of course, the 
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problem is simplified since Ne is then a constant. For horizontal tails 
only partially submerged in the slipstream it is again difficult to assess 
just what portion of the tail is actually affected by the slipstream. 

| In spite of these errors, however, this flight test technique produces 
experimental values of C,, 3 that are in good agreement with theory and. 


should prove to be a valuable tool in determining this parameter. 


III. THE DAMPING IN PTICH, Cy, 


1. Theory: 
Consider an airplane flying along a straight path in equilibrium at 
some particular lift coefficient, i.e., unaccelerated flight. The balance 


of pitching moments is then given by: 


= OV , 
Ome g = Cr, + (=) AC, + Cms Ade (12) 


CVTEK 
At each trim point the lift mst be equal to the weight and hence 


oe is a constant. If the airplane's speed is decreased its lift coeffi- 
cient must increase and the change in pitching moments due to this increase 
must be balanced by the elevator. If the airplane is placed in a steady 
turn at some bank angle the lift must increase in order to provide a vertical 
component of lift equal to the weight. In a turn, however, the lift coef- 
ficient may be thus increased at no decrease in forward speed and the balance 
of pitching moments will be different from those obtained when the lift 
coefficient is increased in level flight. That is, the change in pitching 
moments, that must be balanced by a change in elevator deflection, will be 


different, and therefore the elevator deflection required to produce a given 








eS PU 
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increment in lift coefficient will not be the same as that deflection 

required to produce the same increment in lift coefficient in level flight. 
Obviously, this difference in elevetor deflection is that increment required ; 
to overcome the damping in pitch of the airplane. 


If the damping derivative is written in the form C where (d) is the 


m 
de 
operator iT (t= airplane's time characteristic = <7 ), the 





balance of pitching moments for accelerated flight is given by: 


Cmeg, SG (dom) acy + Cm 4°48 + Cm: MSe (13) 
bts 2 
As shown in the section on elevator power, the elevator deflection 


must balance out any moments acting on the airplane if equilibrium is to be 
obtained. Having previously determined the elevator power, C,, 5? from the 
leg trim curves, equation (13) offers a method of determining Cnig from 
flight test data. a 

The airplane lift coefficient is defined as Crs ae at a 
given velocity, V, then, the lift coefficient on a 2-g curve is exactly 
twice that on a l-g curve. If curves of , versus V,,, for several values 
of g are taken in steady turns and replotted as 6, versus Cr? the increment 
in §, between two of these curves (between l-g and 2-g for example) will be 
that increment in elevator deflection required to overcome the damping in 
pitch of the airplane. If no damping moments existed the curves of § 9 versus 
ke for various g's would coincide since a given elevator deflection would 
produce a given lift coefficient at any acceleration. This assumes, of course, 
that the airplane's static stability does not vary with the different slip- 
stream characteristics encountered at the various power settings used in ob- 


taining the data. 
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If this assumption is made, these curves of 5. versus C, for 
a 
various g's may be reduced to obtain the damping in pitch derivative, 


Ona in the following manner. The airplane's angular velocity in pitch 
may be shown to be: 


6 = a(n - ) (14) 
V =a}: 
(Reference 1.) 


Or, in terms of the operator (d): 


de = fta (1-4) (15) 





For any value of the airplane lift coefficient, then, d@ may be easily 


determined. The equilibrium condition is given by: 


A 


Therefore : 
Cms: AS 
Cm do = = mee eer (47) 


2 Determination of Omg: 


Equation (17) developed above indicates that C can be computed 


Ne 

from flight test data obtained in accelerated flight. The technique may 

be outlined as follows: / 

1. Curves of § @ versus Cr. for seveyel values of n, including n = 1.0, 
are obtained for one or more center of gravity locations of the 
airplane. 

2 For a particular Oi and a particular n, compute d@, the airplane's 
rate of pitch in terms of the operator (d). 

3. At this particular Cr, determine the AS m required to accelerate 


the airplane from n = 1.0 to the particular n under consideration. 








i 


4. Having previously determined Cn s as Sales in the section on 
elevator power, the damping in pitch, Cmge, is then easily de- 
termined from equation (17). 

5. This process is repeated for a number of Cy 's and n's. Although, 


in general, C " is not a function of lift coefficient, it is con- 


d 
venient to plot Omae versus ve in order to establish an average 
value for the derivative. The scatter of the computed points is 
usually considerable. However, if a sufficient number of points 
are taken, the trend of the curve will be indicated well enough 
thet the curve itself can be determined with fairly good accuracy. 


The final value for C then, is the average value of this curve. 


Mae? 
3. Results: 

The curves of § , versus V, obtained in accelerated flight were replotted 
as curves of oF versus om as shown in Figures 6, 7, and 8. From these 
curves Cae was computed for every possible combination of Che and n as oute- 
lined in the preceding section. A plot was then made of Cmag versus Cle as 
shown in Figure 9. As expected, considerable scatter of the individual 
points wes evident in this plot. However, the trend of the curve is clearly 
indicated. While the individual points computed are not shown in Figure 9, 
the range of the scatter is indicated by the cross-hatched band. As Cma6 
is,in general, not a function of lift coefficient, the experimental curve 
of Figure 9 was considered to be codaeene over the range of lift coefficient 


investigated. The average magnitude was determined to be = 0.195. 


AERONAUTICAL ENGINEERING LABORATORY REPORT 160, 


; 











AE 


oe a 





4. Theoretical Calculation of Cnag 

The total damping in pitch of an airplane is the sum of the damping 
contributions of the verious airplane components -- fuselage, wing, hori- 
zontal tail, and the propeller. In the normal configured airplane the 
damping due to the tail is by far the largest factor involved. It is the 
usual practice to evaluate the damping in pitch due to the tail and then 
increase the result by some arbitrary factor to take care of all the other 
contributions to the total damping. This arbitrary multiplying factor is 
normally taken equal to 1.10. 

The damping due to the tail occurs as a direct result of the ehange 


in the effective angle of attack of the tail produced by the pitching 


velocity of the airplane, i.e., the angular velocity 6. This change in oe, 


is given by: ; 
V 


The pitching moment produced by this angle of attack is: 


Cy = — One Vee b 
bil Bey 


Then: 
Cm = -— % We V de 


d(g¢) V 


Dividing each side by the time characteristic, T : 


dCm =, V7 
ae ee te 
From which: d a . 


(18) 


(19) - 


(20) 


(2.1) 


(22) 
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To account for the rest of the airplane's damping in pitch multiply 


by 1.10 giving as the final expression: 


Cmde = — (210 a4 V ne et (2.3) 


Calculation of the theoretical value of the damping in pitch according 


to this equation is as follows: 


1. a = 0.048 per degree = 2.75 per radian 
Ze Vv = 0.428 








3- Assume NH, = 0.90 | | 
Met «6 (19812 8G 

Cc 59.02 a 
De = m 


‘_ pce | 
= © x 12 
= ee5 x 159.20 x 59.02 x 3202 
ae 
'f 


28.0 


1 
B 28.0 — 0.0357 


=~ 1.1 x 2675 x 0.428 x 0.90 x 2.8 x 0.0357 


(©) 
5 


a 0.1165 


5. Conparison of Theoretical and Experimental Values for Cmag: 
The theoretically determined value for the damping in pitch (Crip - 
- 0.1165) is plotted in Figure 9 to facilitate comparison with the experi- 


mental curve. The value of this parameter as deduced from the flight test 
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data is - 0.195 as shown in Figure 9. The experimental value is 66.8% 
higher than theoretical. This large discrepancy between experimental and 
theoretical results may be due to one or more of several factors. In the 
following paragraphs these possibilities are discussed in some detail. 


In computing the experimental values of Oma at the various CL, 's ; 


6 
the magnitude of the elevator power, C,, S$ used in equation 17 to determine 
Ong at a particular Gy. ras taken directly from the experimental curve of 
Figure 5. As pointed out in the discussion on elevator power, the experi- 
mental curve of C,, s versus CL, was also somewhat higher than theoretical.: 


This discrepancy was attributed to the fact that the variation in Ny was 





unknown and to the assumption that CLe = Choe It has been shown that if 





this assumption is not made, that is, the elevator power is computed on the 
basis of the lift coefficient of the wing alone when CL. is not equal to 
Cr» then the elevator power will be reduced to a lower value. (Reference 3) 
Since Omag varies directly with Cos such a reduction obviously would also 
decrease the magnitude of Cy 4ae° 
Another important factor that has been assumed negligible in these 
calculations is the varietion of the airplane's static stick-fixed longi- 
tudinal stability with the variation in power. It was assumed that the 
change in elevator angle required to accelerate the airplane to ann greater 


than 1.0 at a particular Cl» i.e., the A& between n curves, was entirely 


due to the damping in pitch of the airplane. It can be shown that: 


d5e _ _ dCm/dcy (24) 


dc, ' Cms 
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That is, the elevator angle required to vary the equilibrium lift 
coefficient varies inversely with the elevator power, C,, 5°? and directly 
with the static stick-fixed longitudinal stability, dCy/dC,. The as- 
sumption, then, that the total AS, between two different n's at a parti- 
cular Cy includes the assumption that the static stability remains a 
constant throughout the test. This is probably a poor assumption in view 
of the fact that in some airplanes there is a wide variation in dC,/dy, 
with variation in power. (This variation is a function of the reletive 
positions of the airplane's thrust line and its center of gravity location 
in the vertical direction). If this is the case, then, a wide variation 
in power as experienced in these flight tests may be expected to introduce 


considerable error in the value obtained for C If the variation of 


Mg ° 
the static stability with power is known, or can be accurately predicted, 
such that the change in 8, required by the static stability change can be 
deducted from the total AS& e required to accelerate at a given Cy, then 
Crag can be determined more accurately and its magnitude will be lower for 
the case where static stability increases as power is increased. 

A third source of error which may account for pert of the discrepancy 
between theoretical and experimental results lies in the theoretical cal- 
culation itself. It will be remembered that the damping in pitch of the 
airplane was consideree to be due primarily to the horizontal tail damping. 
fal other damping contributions due to the fuselage, wing, and propeller 
were lumped together as a mitiplying factor of 1.10 to be applied to the 
horizontal tail damping. There is some evidence to indicate that this 


factor of 1.10 may be too low for airplanes having low relative densities. 
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Zimmerman developed a theoretical formula for Cmae in 1935 in which this 
multiplier was taken to be 1.25. (Reference 4.) Since this time the 
relative density of most airplanes has increased considerably and the factor 
of 1.10 is generally considered adequate. For modern light airplanes with 
low wing loadings, however, the relative density ee » remains low, 
which in itself tends to increase the damping parameter, Cmag, since Cmae 
veries cameeeey" Fao Nevertheless, it is believed that increasing 

the mitiplier from 1.10 to 1.25 or higher when working with light airplanes 
may lead to more accurate results in the theoretical calculation of Crag: 
Increasing the mltiplying factor to 1.25 yields a Cmag for the Cessna 140 
of ~ 0.133. The experimental value obtained from the flight tests is then 
only 46.6% higher than theoretical as compared with 66.8% when the factor 

is 1.10. This new value, of course, is calculated under the same assumptions 
as before, i.e., the variation in C, 5 and the change in static stability 


due to power are neglected. 


€. Evaluation of Experimental Technique for Obtaining Cmasp: 
In the preceding sections a flight test technique for evaluating the 


airplane's damping in pitch has been outlined and discussed in some detail. 
The usual question now arises as to whether or not the method produces 
results of sufficient accuracy to be practical. 

The pilot Lebtinigte required to obtain trim curves in accelerated flight 
is relatively simple. The steady banked turn at constant speed is the 
easiest and most precise method of securing this data. The variables involved 
-- the airspeed, acceleration, and the elevator deflection -- can be measured 


with considerable accuracy. If the elevator power is accurately known, then, 
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ae should be possible to reduce the flight test data to a reasonably 
exact value for Omi9° 
The major problems to be solved in calculating Cnae in this manner are 
the variation of Cn with lift coefficient and the variation of static 
stability with power. The sources of error in determining os 5 have already 
been analyzed in previous sections. The static stability variation is parti- 
cularly difficult with propeller driven airplenes. If the aircraft is jet 
or rocket powered there may still be i~portant variations in the stability, 


but usually they are far less serious then for propeller driven airplanes. 


In high speed aircraft, however, there are additional problems arising from 





inaccuracies in the measurement of the extremely small elevator deflections 
required. 

In spite of the difficulties encountered in separating the variables and 
minimizing the errors encountered, the value of this technique for obtaining 
the damping in pitch derivative should not be overlooked. In the past it has 
been necessary to rely on avproximations of this parameter. At best, these 
theoretical approximations were in the nature of an "educated guess". Recently, 
more elaborate techniques for obtaining the longitudinal derivatives iT 


frequency response of the airplane have been developed. OC can be determined . 


Mae 
in this manner but the method requires precision instrumentation and data 
reduction. Errors are cumulative throughout the procedure such that by the 
time the value of the damping derivative is reached the results are likely to 
be much less accurate than those obtained from the steady state flight tests. 
It is believed, therefore, that with good pilot technique in obtaining the re- 


quired data, dependable instrumentation, and the elimination of as many of the 





AERON 


7 
AUTICAL ENGINEERING LABORATORY REPORT 160 


: 


——e a 


variables as possible, the steady state flight tests for determining the 
damping in pitch derivative as outlined will prove to be the most accurate 
method of obtaining this derivative for subsonic aircraft. Some refinements, 
such as a means of handling the static stability changes with power, are 


desirable and should be the subjects of future studies. 


CONCLUSIONS 
From the results of the flight test program conducted on the Cessna 140 
airplane, it is concluded that: 
1. The elevator power, Cp, §& » can be predicted with good accuracy from the 


unaccelerated trim curves for the airplane provided that: 





a). Power effects are carefully analyzed in order that errors due to 
variation in the tail efficiency, Nt» can be minimized. 

b). Consideration is given to the fact that the airplane lift coefficient 
may not be zero when the wing lift coefficient goes to zero. The as- 
sumption that Che = Cr is fairly good in most cases but the accuracy | 
of the data reduction will be increased somewhat if the assumption 
is not made and the calculations are corrected accordingly. 


2 The damping in pitch, © g? can be predicted with good accuracy from the | 


ng 
accelerated trim curves for the airplane provided that: my ¥ 
a). The elevator power, Cn g§ » can be determined accurately. 
b). The variation in stetic stick-fixed longitudinal stability with power 

is either small enough to be neglected or can be accurately predicted. 
the miltiplying factor of 1.10 


wo (=) 
used to correct the horizontal tail damping to include the damping of the 


3. In making the theoretical calculation of © 


whole airplane may be too low for airplanes having low relative densities. 
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REC OMMEN DAT IONS 


On the basis ef the above conclusions the following recommendations 
made : 


The unaccelérated trim curves of the Cessna 140 should be further investi- 


gated. In particular, the l-g curve for the forward center of gravity 


location should be re-flown to check it below cr. - 0.4. The behavior 

of this curve should be carefully studied. 

A further study of the elevator power of the Cessna 140 should be made 

in which C, § is computed on the basis of the wing lift coefficient alone 
as outlined in Péfereiieer 3. 

A flight test program should be carried out on the Cessna 140 to determine 
the yariation of static stick-fixed longitudinal stability with power in 


order that the computed value of CG can be corrected for this effect. 


Hide 
A further study should be made in the theoretical determination of Cna0 
to investigate the effect of low relative density of the airplane on the 

result. Damping in pitch due to the propeller, wing, and fuselage should 
be studied to determine more accurately what portion of the total damping 


must be attributed to the components. 
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APPENDIX I 


LIST OF SYMBOLS AND ABBREVIATIONS 


See eee eee 


fir density, slugs/cu. ft. 


True airspeed 


= 


Indicsted eirspeed 


<j 


Calibrated airspeed 





Cc 
Ve Ecuivalent airspeed 
a Elevator deflection, degrees 
9 Dynamic pressure, lbs/sq.ft. = hove 
Fe Acceleration of gravity, 32.2 ft/sec/sec. 
n Normel acceleration in units of g. 
Cc Mean aerodynamic chord 
S) hrea 
W Weight 
m Mess = W/g 
Mi Moment 
L Lift 
Cte Airplane lift coefficient 
CL Wing lift coefficient 
By, Slope of horizontal tail lift curve 
Oy Slope of wing i1ift curve 
X. Angle of attack of wing 
Ky Angle of attack of tail 
= See” 
V Tail volume coefficient = ie Sian 
he Tail length (airplane c.g. to aerodynamic center of tail 


N+ Tail efficiency factor = slipstream q/ free stream q 








@ 


QR. oe oS gc 
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Blevator effectiveness dX, /d §, 





Airplane's time characteristic x a 
Airplane's relative density »s i a 
Angle of pitch 


Anguler rate of pitch 


Operator = OCA. 
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APPENDIX JIT 


WEIGHT AND BALANCE COMPUTATION 





Seale Calibration 


Scale Load in Pounds 
fs. 3 =e = [= [= =e] 
pan [sniffs] | o[ olin 
pam [aspespolerl = [= Dale | peopel 
Fas [ool ow [aw [om [ow [om wales] | mm fe 


Weight* and Balance 


Condition : Oi1 and Fuel Tanks Full 


Test Equipment Installed 


1. Weight: 












Seale Weight 
co 


Right Wheel oi! 
on wef D| : 






* Airplane weighed in level flight attitude. 
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2e Balance: 
Datum line (reference) ......+ +--+. Root leading edge 
Mean aerodynamic chord . . . . « « « « « « « © « « 5902" 


Root leading edge to mean aerodynamic chord 
leading edge. . .. 0.47" 


Datum line to main gear contact. . . 6 +» « e« « «© « 3631" 


Datum line to tail wheel contact .......-. 201.10" 


a) Normal center of gravity location: 


=, (Wp,+ Vp) Dy + We Dy (527 + 499) (3.3L) + (74) (201.1) 
Gross weight 1100 


x' ss 16.64" 


x! Dist. root L.E. to ma.c. LeE. 2 16.64" = 0.47 = 16.17" 


m | 
In 


x - 50902) s 27 he MeGBeCe 


b) Weight of two pilots and equipment = 335 lbs. 
Located at normal center of gravity location. 


Weight of airplane and crew = 1100 + 335 s 1435 lbs. 


c) Forward center of gravity location: 


Added weight: 16 lead weights = 57.0 lbs. 26.5" forward of datum line. 


a 14,92 oe 


x ad 14.99 = O47 = O. i ad 6% eadeCe 
TS 59.02 ae O8 Beets 
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ad) Mid center of gravity location: 





Added weight: 1 lead weight = 1.5 lbs., 191" aft of datum line. 


°S - (1435 >.< 16.64) + (5 x 191.0) = 16.80" 
1436.5 

-_ O pe / tt 

x, = 16 sar | - 27.6% MeOBcCe 


_e) Aft center of gravity location: 
Added weight: 4 lead weights = 15.5 lbs., 191" aft of datum line. 


—- 6,6 15.5 x 191,0 . 
X3 = 1,50 = 18.558 


x = 18. ee vi =~ edeCe 
x, = ese el - 30.6% mea.c 
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